Two-point phase correlations of a one-dimensional bosonic Josephson junction 
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We realize a one-dimensional Josephson junction using quantum degenerate Bose gases in a tun- 
able double well potential on an atom chip. Matter wave interferometry gives direct access to the 
relative phase field, which reflects the interplay of thermally driven fluctuations and phase locking 
due to tunneling. The thermal equilibrium state is characterized by probing the full statistical dis- 
tribution function of the two-point phase correlation. Comparison to a stochastic model allows to 
measure the coupling strength and temperature and hence a full characterization of the system. 



Josephson dynamics between weakly coupled macro- 
scopic wave functions have been observed in supercon- 
ductors [HH], superfluid Helium [31 0], an d recently using 
Bose-Einstein condensates in double well potentials [5]- 
[7]. The bosonic Josephson junction (BJJ) is especially 
interesting, as particle interactions lead to additional dy- 
namical modes such as quantum self trapping or tt phase 
modes [31 [5] and finite temperature leads to enhanced 
fluctuations of the observables [S]. In contrast to other 
implementations, the BJJ enables complete experimental 
control over all relevant system parameters such as the 
coupling strength or relative population together with 
direct access to the conjugate observables number and 
phase. Theoretical work has mostly employed a two- 
mode approach to describe the finite temperature equi- 
librium system and dynamical properties 110) . 

One-dimensional (ID) Josephson junctions show a sig- 
nificantly enriched physical behavior, as the two involved 
wave functions can not be described by single quantum 
modes any more. The non-interacting ID junction rep- 
resents an implementation of the Sine-Gordon Hamilto- 
nian which occurs in widespread areas of physics [111 I12j . 
In the ID bosonic junction interactions and finite tem- 
perature are expected to cause dynamical instabilities of 
the classical Josephson modes [T3J . Whether quasi-static 
phenomena such as quantum self-trapping persist in ID 
is issue of ongoing discussion [T3] . 

In this work we realize and fully characterize a one- 
dimensional bosonic Josephson junction using quantum 
degenerate Bose gases in a tunable double well potential. 
The finite temperature equilibrium state is marked by 
the competing effects of thermally driven phase fluctua- 
tions and phase locking due to tunnel coupling. Fluctua- 
tions of the relative population are < 1 % and can be ne- 
glected [5] . We probe the coherence properties of the cou- 
pled system by performing matter wave interferometry. 
Comparing the statistical distribution function of two- 
point phase correlations to a stochastic model [10l [15] , 
we measure the coupling energy or the temperature of 




FIG. 1. Two coupled Bose gases are released from a dou- 
ble well potential. The matter wave interference pattern 
emerging during expansion gives access to the relative phase 
tp{z) along the samples. Fluctuations in the absolute phase 
of the two Bose gases transform into visible density fluctua- 
tions [161 117] . We characterize two-point phase correlations 
of the system by measuring the statistical properties of the 
difference of relative phases Aip(z) = ip(z) — tp(z'). 



the system. 

The experiments are performed in a horizontally ori- 
entated double well potential, generated on an atom chip 
using radio-frequency (RF) induced adiabatic states [HI 
[19] . Different double well parameters and hence dif- 
ferent coupling strengths are realized by using differ- 
ent RF amplitudes, with distances between the minima 
ranging from 1.2 fmx to 2 /xm and a barrier height be- 
tween h x 1.5 kHz an h x 12 kHz. The trap frequencies 
of the individual wells are measured as v y — 2.9 kHz 
and v x = 3.3 kHz along the tightly confined directions 
and v z = 18 Hz along the longitudinal direction. The 
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two samples are coupled along the strongly confining di- 
rection y. Note that the double well potentials remains 
static throughout an experimental sequence. 

Starting from a thermal gas of 87 Rb atoms in the 
\F = l,nip = —1) state, we create a system of coupled 
degenerate Bose gases of adjustable temperature using 
forced evaporative cooling in the static double well po- 
tential. This, together with a phase of plain evapo- 
ration of at least 180 ms, ensures a system in thermal 
equilibrium. Each well contains 3200(510) atoms, corre- 
sponding to an in-situ peak line density of 72(7) ^m _1 , 
and a chemical potential of fi/h — 1.72(0.17) kHz [20] . 
With temperatures of typically 150 nK, a regime where 
p < UbT < hh> XtV is realized. As the physics consid- 
ered here concerns only long-wavelength excitations, each 
sample can be regarded as one-dimensional [501 HI] E3 • 

To probe phase correlations of the BJJ we perform 
matter wave interferometry. After suddenly (< 10 /zs) 
switching off all confining potentials the two samples 
overlap in 46 ms time-of-flight (TOF) expansion. The 
emerging interference pattern is recorded using a light 
sheet fluorescence imaging [23] . The relative phase 
(p(z) = f L {z) — (fi R {z) between the left (L) and right (R) 
Bose gas is determined by fitting a cosine function with 
a gaussian envelope to the density profile in each of the 
4/xm wide pixel slices (see Fig. [I]). The fitting introduces 
an uncertainty 5<p(z) on the value of (p(z) which follows 
a normal distribution with zero mean and a standard de- 
viation ranging from 0.1 to 0.5 radians. The analysis is 
restricted to a maximum distance of 40 pm over which 
the atomic density drops by 30%. The experiment is re- 
peated 500 times for each parameter set to provide a full 
analysis of the statistical properties of f(z). 

A single trapped ID Bose gas at finite temperature is 
characterized by a density profile po(z) and a spatially 
fluctuating phase [H] responsible for the exponential de- 
cay of the first order correlation function on the length 
scale At = 2ft 2 'po/ '(m^T) [31]. Interference of two in- 
dependent samples can be used to probe the coherence 
properties of single ID Bose gases [25Tf27] . 

Introducing a tunnel coupling strength J between 
two ID Bose gases implements a combined system with 
an additional degree of freedom in the relative phase 
(p(z) [10] [28]. It is associated with the additional length 
scale I j = y/h/AmJ, which represents the typical dis- 
tance on which tunnel coupling restores a spatially con- 
stant relative phase. For lj < At the system becomes 
phase-locked, i.e. the coupling counterbalances the ran- 
domization due to thermal excitations. 

Two-point phase correlations along the coupled ID 
samples can be probed through the statistical proper- 
ties of the difference of the relative phase Aip(z) — 
tp(z) — (p(z') where z = \z — z'\. The scaling of the distri- 
bution of Aip(z) with z gives access to the spatial exten- 
sion of the relative phase correlations and hence to the 
relevant length scales At and lj [10] . A narrow distribu- 
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FIG. 2. Distributions of the difference of the relative phase 
Aip(z) for decreasing tunnel coupling (from left to right, com- 
pare table [I]). (a,b) High coupling yields a narrow distribution 
and high coherence over the entire length of the system, (c-e) 
Decreasing coupling leads to an increasingly fast loss of spa- 
tial phase correlations and thus a randomization of relative 
phases. A slight bias of the distribution can be assigned to 
a small relative velocity (20/im/s) of the two samples during 
the expansion. 



tion of Aip(z) (peaked around 0) indicates high coherence 
and phase locking, whereas a broad distribution (between 
— 7r and 7r) characterizes uncorrelated phases. Figure [2] 
depicts measured distributions of A(p(z) for different dis- 
tances and tunnel couplings, ranging from strongly phase 
locked (a) to almost independent samples (e). 

Calculating the real part of the phase correlation func- 
tion C(z) = Re (e^-^'A = (cos(A^(z))) allows to 

quantify the spread of the distribution of Aip(z) [15]. By 
changing from 1 to 0, C(z) describes the transition from 
spatially locked to uncorrelated relative phases for in- 
creasing z. The characteristic length scale and shape of 
this decay is determined by the BJJ parameters po, J, 
and T. Figures [3] and [4] show measured values of C(z) 
for different values of J and T at constant atom number. 

To compare our results to theoretical predictions of [TO] 
we simulate single realizations of tp(z) using an Ornstein- 
Uhlenbeck (OU) stochastic process [15] . These real- 
izations exactly follow the phase correlations predicted 
in |10) for uniform systems but allow us to account for 
a finite number of experimental realizations and finite 
system size along with a spatially varying density po(z) 
(calculated according to [3D]) and a density-dependent 
coupling J(po(z)) within a local density approach [2"5] . 
From this we construct single realizations of the in-situ 
wave functions i/) ' R (z) = \J po(z) exp ( iipg ' (z) ) each 
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FIG. 3. Real part of the phase correlation function C(z) for 
different couplings and T ~ 155 nK (compare Table [B. The 
symbols represent experimental values, derived from the data 
shown in Fig. [2] A stochastic OU model allows to estimate 
the couplings J\ to Js (compare Table [T| . The colored areas 
depict two standard deviations on the values of C(z) obtained 
by repeating the analysis 30 times. A simulation of C(z) 
for completely uncoupled Bose gases is displayed as a black 
dashed line and indicates the limits of our imaging resolution. 



representing single outcomes of the experiment. 

When comparing to experimental data, we have to in- 
clude the effects of the expansion and the detection pro- 
cess. Since the expansion of the system can be consid- 
ered ballistic [TH1 [T7] we numerically compute \tfj L ' R ^ — 



U(t) 



V-'o 



for a TOF t = 46 ms, where U(t) is the 



free evolution operator. The resulting relative phase 
ip(z) — ip L (z) — tp R {z) and the density p(z) after ex- 
pansion are then individually convoluted with the point- 
spread-function of the imaging system, which we estimate 
to be gaussian with an RMS width of 7 pm [23] . Finally 
the uncertainty Scp(z) introduced by the fitting procedure 
is added to the resulting relative phase, coarse grained to 
the 4 /jm pixel size. 

Averaging over 500 realizations of the process as in the 
experiment, the experimental results can be reproduced 
by adjusting the simulation parameter po, Jou, and 7b u- 
Repeating the described procedure many times allows to 
evaluate the uncertainty on the estimation of C(z) due 
to the finite number of experimental realizations. 

Figure [3] compares measured and simulated values of 
C(z) for five different tunnel couplings (J\ to J5), with 
fixed temperature Ttof and density profile po(z). The 
temperature Ttof is determined by fitting a Bose func- 
tion to the thermal wings of the recorded density pro- 
files [23]. To match the OU simulations of C(z) with 
the experimental data, only the parameters J\ to J5 are 
adjusted (see Fig. [3] and Table [jy). 

This correlation analysis provides a direct experimen- 
tal measure of the coupling strength in a BJJ, without the 
necessity to probe dynamical properties or rely on sim- 
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FIG. 4. Real part of the phase correlation function C(z) for 
fixed coupling J3 and different temperatures Tj to T3. Com- 
parison with the stochastic simulation allows thermometry of 
the ID bosonic Josephson junction (compare Table |l|. 



ulations of the system in the double well potential. We 
show here that our experimental control allows to adjust 
the tunnel coupling over two orders of magnitudes. 

A comparison with an alternative numerical simulation 
based on the time-dependent ID Gross-Pitaevskii equa- 
tion shows, within a factor of two, a good agreement with 
our results (compare Table [I] and see J£j ) . Note that a 
factor of two in tunnel coupling corresponds to less than 
3 % uncertainty in parameters defining the double well 
potential, highlighting the need for a direct and precise 
experimental measurement of </, as realized in this work. 

In cases where the tunnel coupling is known, the anal- 
ysis of C(z) can be used for thermometry of the Joseph- 
son junction as shown in Fig. [4] [9]. Here, the trapping 
potential remains unchanged and J is fixed to J3. The 
temperature T is set by changing the final position of 
the cooling RF-field while keeping the total atom number 
constant. The temperatures Tou obtained by adjusting 
the OU simulations to the data are in excellent agree- 
ment with the independently measured Ttof (compare 
Table [I]). For experimental systems featuring an optical 
resolution better than the thermal coherence length, the 
analysis presented in this work would even allow the si- 
multaneous determination of the coupling J and the tem- 
perature T. 

The experimental data as well as the OU simulation 
give access not only to expectation values but also to 
the full distribution of A(p(z). Figure [5] shows exemplary 
distributions of Aip(z) and its cosine, for z = 32 pm and 
couplings Ji, J3, and J5. The excellent agreement be- 
tween experiment and simulation indicates that we have 
indeed realized an equilibrium system as decribed in [lOj . 

In summary we have used two-point phase correlations 
to fully characterize a tunable ID bosonic Josephson 
junction in thermal equilibrium. Comparing the experi- 
mental data to the results of a stochastic model based 



4 



Fie- 3 


TTrwfnTO 
- 1 1 Ur ^u.xv J 






^ \r i,,t 'J 




j 

J\ 


104^0 J 


ID. 0-^1 




O.o 


O O Q Q 


Jl 


150(5) 


5-7.5 


4.8 


5.4 


4.9-6 


Js 


153(5) 


0.65-0.9 


0.15 


5.3 


14.4-16.8 


J 4 


154(5) 


0.05-0.15 


0.05 


5.1 


34.9-60.4 




163(5) 


<0.1 


<0.05 


4.8 


>47.8 


Fig. 4 


Ttof(iiK) 


Tou(nK) 


Mm 




lj(pm) 


Ti 


155(10) 


125-180 
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TABLE I. Comparison of experimentally measured and sim- 
ulated parameters. All spatially dependent values are given 
at the position of the peak atomic density (z = 0). 
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FIG. 5. Distributions of Atp(z) (a-c) and cos(Aip(z)) (d-f) 
for z = 32 fim and three different couplings (Ji, J3, J5). Red 
histograms represent measured, blue lines simulated distribu- 
tions for 500 realizations, (a-c) The mean value of Atp{z) 
(dashed gray line) does not change while the phase spreads, 
(d-f) The mean value of cos Aip(z) (dashed gray line) decays 
from about 1 (strongly coupled) to (uncoupled). 



on an Ornstein-Uhlcnbeck process allows us to deter- 
mine the coupling strength or the temperature of the 
system. This full characterization marks the starting 
point for further studies of dynamical properties and non- 
equilibrium states of coupled one-dimensional systems. 
Further research will investigate the interplay of spatial 



density and phase correlations. 
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